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Abstract. In this study, an adaptively hybrid method was proposed to improve the performance of fractal coding 
methods. First, we found that the range blocks with large variances (RBLVs) play a crucial role in degrading 
decoded images, and the effect of the remaining range blocks with small variances (RBSVs) can be ignored. 
Second, RBLVs were designed to be encoded in an extended domain block pool (EDBP), and the remaining 
RBSVs were encoded with the no-search fractal encoding method. Moreover, an effective method to adaptively 
divide the range blocks into the above two categories was proposed. Finally, four fractal coding methods were 
adopted to assess the performance of the proposed method. Experimental results show that, compared with the 
previous methods, the proposed method can achieve better-decoded image quality with fewer bits per pixel and 
fewer computations. 
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1 = Introduction 


Different from conventional image compression techniques, fractal image coding is 
characterized by its novel idea, potential high compression ratio, fast decoding, and 
resolution independence.!'“! Thus, it has attracted much attention from researchers globally. 
However, owing to its high computational complexity in the encoding process, fast fractal 
image coding becomes one of the most important topics in this field. In recent years, many 
researchers have attempted to overcome this problem by converting exhaustive searching into 
local matching.>-”! In this case, as the number of block-matching operations decreases 
significantly, the encoding process can be done in a shorter time. To further accelerate the 


encoding process, no-search fractal image encoding should not perform the block matching 
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operations, and a domain block is directly appointed as the best-matched domain block. Thus, 
no-search fractal image encoding can achieve real-time encoding and a higher compression 
ratio at the expense of the quality of the decoded images."! In addition to image compression, 


fractal coding methods have also been gradually implemented in other image processing 


[9-15] 16-18] 19,20] 


applications such as image denoising, image retrieval,| image hashing,| 
watermarking,’ 1-23] image super resolution!*+7*!, and MRI images.!??! 

In this study, we found that the range blocks with large variances (RBLVs) play a 
crucial role in determining the quality of decoded images. However, the effect of the 
remaining range blocks with small variances (RBSVs) can be ignored. Therefore, we 
designed an extended domain block pool (EDBP) for RBLVs. Because EDBP can effectively 
reduce collage errors, the quality of decoded images can be expected to improve. 
Furthermore, to accelerate the encoding process and reduce the bits per pixel, we adopted the 
no-search fractal coding method to encode RBSVs. Moreover, the bottom limit of the actual 
percentage of an accumulated collage error (BLAPACE) was proposed to adaptively divide 
the range blocks into the above two categories. Finally, four fractal encoding methods, 
Jacquin’s,"! Chaurasia’s,! Zheng’s!®!, and Gupta’s methods,!7! were adopted to assess the 
effectiveness of the proposed method. Experimental results show that compared with the 
previous methods, the proposed method can achieve a better quality of decoded images with 
fewer bits per pixel and fewer computations. 

The remainder of the study is organized as follows: an introduction about conventional 


fractal image coding is reviewed in Sec.2. In Sec.3, a hybrid fractal encoding method is 


proposed. By analyzing the importance of variances, RBLVs are encoded in EDBP, and 
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RBSVs are encoded with no-search fractal image coding. An adaptive method is also 
proposed to divide the range blocks into the above two categories. In Sec.4, the experimental 
procedures and performance of the proposed method are described and discussed in detail. 


Finally, the conclusion is given in Sec.5. 


2 Conventional fractal image coding 


Fractal image coding aims to establish an iterated function system (IFS) whose fixed point 
can well approximate the input image. It mainly consists of a series of block-matching 
operations between range and domain blocks. First, the range blocks divide the MXN input 
image uniformly into nonoverlapping BxB blocks; Ri, i=1,2,3,..., NumR, and domain blocks; 
D,, j=1,2,3,..., NumD, can be obtained by sliding a 2Bx2B window over the input image from 
left to right and top to bottom, where NumR and NumD represent the total numbers of range 
and domain blocks, respectively. In this study, the sliding step 6 was set as 2B. To make 
domain blocks match range blocks well, domain blocks are first contracted to the same size 
as range blocks and then extended with eight isometric transformations, which consist of the 
identity, rotations through 90°, 180°, and 270°, and reflections about the vertical axis, 
horizontal axis, first diagonal, and second diagonal. Figure 1 illustrates the eight 
transformations of an image. After performing affine transformations on the contracted 
domain blocks for arbitrary range blocks, the corresponding best-matched domain block can 


be obtained by minimizing the following function 
CE(R)= nmin [min [lk -sD arf) l 
D 5,0 B 


where J denotes a B x B matrix whose elements are all ones. s and o denote the scaling and 
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offset coefficients for the affine transformation, respectively. CE(R) is the collage error of the 
range block R, and D denotes the domain block in the conventional domain block pool 
(CDBP). With the least square method, we can calculate s and o as 
= =H. = 
s=(R-71, D-dt\/|p-dil , o=r-sd 2 
where (*.*) denotes the inner product. 7 and d are the mean values of R and D, 


respectively. Figure 2 illustrates the above fractal encoding process. 


At the decoding phase, an arbitrary image with the same size as the input image can be 
selected as the initial image. The same transformations from the best-matched domain to the 
corresponding range blocks as those in the encoding process were applied recursively to the 
initial image. After approximately 10 iterations, the decoding process converged into the final 
decoded image. As an example, we used 256 x 256 Peppers as the input image. If the rotated 
256 x 256 Peppers were chosen as the initial image, as shown in Fig.3(a), the first five 


iteration images obtained in the decoding process are shown in Fig.3(b)-(f), respectively. 
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Fig.1. Eight isometric transformations. (a)-(h) represents the identity, rotations through 90° 180°, 270°, and 


reflections about the horizontal axis, vertical axis, first diagonal, and second diagonal, respectively. 
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Fig.2. Illustration of the fractal encoding process. 
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Fig.3. Fractal decoding process. (a) Initial image. (b)-(f) First five iteration images in the decoding process. 


3 Proposed method 


The peak signal-to-noise ratio (PSNR) was adopted in this study to measure the quality of 


decoded images as follows: 


PSNR = log, (ss a sl a 4 n= P (i, D 3 


P 


where fO™2'™! and f°eod4 denote the original and decoded images, respectively, and M and N 
are the height and width of the images, respectively. Moreover, we defined the average 
collage error (ACER) as the average of the collage errors of all range blocks, which can be 


formulated as 


NumR 
ACER = >) CE(R,) Name á 


i= 
where NumR denotes the total number of range blocks. Furthermore, a logarithmic 
relationship between ACER and PSNR of the decoded image was established as? 

PSNR = £, + B,log,, (ACER) 3 
where fı and f2 are constant coefficients. Because f2 is a negative value, we know that a 
smaller ACER implies better-decoded image quality. Eqs.* 5 imply that the PSNR of the 
decoded image can be directly determined by the sum of all collage errors. 


For arbitrary R, by substituting ? into ', we have 
1 
CE(R) = -7R -sD - ol | 
1 2: 2 — 1/2 
==, (|R-Af -sp-a f) 


[e-m.0-ay 


1 2 
Seal 
B? Rez |R -F1 |p-ar| 


= SelR-7f (I-Lcc*)< e-a f 
where LCC denotes the linear correlation coefficient between R and D. From Eq.°, we can 
see that for an arbitrary R, its variance provides a top limit for its collage error, and the range 
blocks with large collage errors should have correspondingly large variances. Thus, only 


RBLVs can lead to large collage errors. If we use LCCmin to denote the minimum of the 


LCCs of all range blocks, we obtain the maximum of CE(R) as 


CE(R) == |R-aAf (1-(LCCy,)) 7 


Max 
Combining Eqs.® ’, we have the following inequality 

CE(R) < CE(R) 8 
In the encoding process, we sorted the range blocks by their variances from largest to 


smallest and encoded each range block separately. Subsequently, we obtained the 


accumulated collage error (ACE) as 
Num! 


ACE= Ý CE(R?"*) 9 


i=l 


where Num1 denotes the number of coded range blocks. From Eqs.®°, we have 


ACE < ACE yax 10 
where 
Num! Num1 j ; 
ACE, = 5 CERO), n, = 5 Jro -7| (1-(Lecge) $ 
ml i=l 


where LCCH! denotes the minimum of the LCCs of the coded range blocks. Similarly, for 


the remaining uncoded range blocks, the corresponding ACE can be defined as 


Num2 
ACE Uncoded = > CE( REE) 12 


i=l 
where Num2 denotes the number of uncoded range blocks. Similar to Eq.!°, we obtain 


Uncoded Uncoded 13 
ACE = ACEi 


where 
Num2 Num2 J , 2 7 2 
Uncoded ___ Uncoded _ Uncoded _ =Uncode es Uncoded 14 
ACE vax a > CE(R, D = > |r r 1| (1 (ECCS ) ) 
i=l i=l 


where LCCH denotes the minimum of the LCCs of uncoded range blocks. Moreover, 


7 


202309.00100v1 


chinaXiv 


we have 
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ECC h LCCN 


Min 


By replacing LCCy""" with LCCmin, Eq. (14) can be rewritten as 
. Num2 Num2 2 > 
ACE = > CE( Roan ek as yy | Roe = Pama I| (1 (LCC un ) ) 16 
i=l i=l 


According to Eq. '° and by comparing Eqs. !4 !6, we have 


Uncoded Estimated_1 17 
ACES ACE 


Max 
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(c) 
Fig.4. Schematic of the LCCs for the three test images. 
Generally, for one input image, there are a large number of range blocks. For example, if 


the size of range blocks is set as 4 x 4, 256 x 256 images consist of 4,096 range blocks. We 


selected three 256 x 256 test images, Barb, Baboon, and Man, for example, and Figures 4(a), 
(b), and (c) illustrate the distributions of the LCCs for the three test images, respectively. 
Clearly, following the encoding order, the LCCs distribute randomly. If sufficient range 
blocks were encoded, the minimum of the LCCs of all range blocks, LCC , can be 


approximated with the minimum of the LCCs of coded range blocks, LCC\" , and we 


obtain 


LGC LCG = 18 


Min 


Coded 


By replacing LCCmin with LCC," , Eq.'® can be rewritten as 


INGE ee 5 CE( RE Woe RUncoded _ Uncoded y i (1 S e )) 19 
i=l 
According to Eq.'8 and by comparing Eqs.!* 1°, we obtain 
ACRES 1 3, Cp Estate 2 A 
According to Eqs.'”*°, we have 
ACR URoded < iNet 21 


If we define the actual percentage of accumulated collage error (APACE) as 


APACE = m oe vy) 
ACE + ACE" 
By Eqs. 7!:”*, we can get the BLAPACE as 
ACE ACE a 
ACE at INCI E 7S ACE ae ACK = 
BLAPACE APACE 


From Eq.”*, we know that in the encoding process, when BLAPACE reaches a large 
percentage, such as 95%, APACE will be more than 95%, and this implies that the coded 


range blocks contribute more than 95% of the total collage errors. If we can reduce the 
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collage errors of the coded range blocks effectively, the total collage errors of all range 


blocks will decrease as well. 
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(c) 
Fig.5. APACE and BLAPACE versus the percentage of coded area for the three test images, respectively. 

Figure 5(a), (b), and (c) show the BLAPACE and APACE versus the percentage of 
encoded area for the three test images in Fig.4. We can easily observe the following key 
points: 

1) Coded RBLVs contribute most of the collage errors. From Fig.5, we can observe that 
in the sorted range block sequences, the first 50% of range blocks contribute 97%, 89%, and 
93% of the total collage errors for the three images. This implies that RBLVs play a more 


crucial role than the remaining RBSVs. Thus, we divided the range blocks into two categories: 
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RBLVs and RBSVs, which contribute most and a little of the collage errors, respectively. 
Furthermore, we adopted two different strategies to encode the above two types of range 
blocks: a) For RBLVs, their corresponding best-matched domain blocks were searched in the 
EDBP, which consisted of the CDBP and supplementary domain block pool (SDBP). Besides 
the CDBP as in Sec.2, we also designed an SDBP, which can be constructed with the 
following two steps: 1. Slide a 4B x 4B window over the input image from left to right and 
top to bottom, and the sliding step was set to 4B. 2. Contract all the domain blocks to the size 
of B x B by averaging every 4 x 4 pixels, and perform eight isometric transformations on 
them. For arbitrary RBLV, its best-matched domain block was searched in EDBP. Compared 
with CDBP in Sec.2, EDBP can provide more domain blocks as the candidates for the 
best-matched domain block, and the corresponding collage error can be expected to decrease 
effectively. b) For RBSVs, we adopted the no-search encoding method to encode them. As in 
Fig.6, the domain block with the same center as the range block was directly appointed as the 
best-matched domain block. From Eq. (6), we know that small variances will provide upper 
limits for their corresponding collage errors. This implies that small collage errors can be 
maintained. However, the no-search method required fewer bits per pixel and computations. 
Finally, for the whole image, compared with the conventional method in Sec.2, we can expect 


to obtain better-quality decoded images with fewer bits per pixel and fewer computations. 


Best matched 
domain block 


Range 
block 


Fig.6. Range block and its associated best-matched domain block. 

2) Consistent with Eq. 7, BLAPACE is often smaller than APACE. When BLAPACE 
reaches a large percentage, such as 95%, APACE should be larger than 95%, and this implies 
that the coded range blocks contribute more than 95% of the total collage errors. The range 
blocks whose BLAPACE reaches 95% of the total collage errors were classified as RBLVs, 
and the remaining ones were considered RBSVs. As shown in Fig.5, when BLAPACE 
reaches 95%, the coded range blocks comprise 56%, 73%, and 70% of the total number of 
range blocks. This implies that for different images, the numbers of RBLVs differ 
significantly for the three images. Thus, by Eq. (23), the range blocks can be adaptively 


divided into the above two categories. 


4 Experiment 


In this Section, five-256 x 256 images (Fig.7, boat, peppers, couple, airplane, and lake) were 
selected. The sizes of range and domain blocks were set at 4 x 4 and 8 x 8, respectively. For 
RBLVs, we used 3, 5, and 7 bits to quantize their isometric, scaling, and offset coefficients, 
respectively. Eleven bits were used to represent its position information. For RBSVs, we used 
2 and 8 bits to quantize the scaling and offset coefficients, respectively. Four fractal coding 
methods, Jacquin’s,'"! Chaurasia’s,°! Zheng’s'*!, and Gupta’s methods,'’! were adopted to 
assess the performance of the proposed method by the quality of decoded images, bits per 
pixel, and the percentage of the required computations. The detailed procedures for the 
proposed method are as follows: 


Step 1: For a fractal coding method, given an input image, divide it uniformly into B x B 
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range blocks. Calculate the variances of all range blocks and sort them by their variances 
from largest to smallest. 

Step 2: Establish CDBP and SDBP as in Sections 2 and 3, respectively, and then EDBP 
can be obtained by combining CDBP and SDBP. 

Step 3: Encode the range blocks separately, and their corresponding best-matched 
domain blocks are obtained by searching within EDBP. By Eq. (23), when BLAPACE 
reaches 95%, the coded range blocks are classified as RBLVs, and the remaining RBSVs are 


encoded with a no-search fractal encoding method. 


(d) Airplane (e) Lake 


Fig.7 Five images used in the experiment. 

Table 1 illustrates the comparison results between Jacquin’s and the proposed methods. 
Clearly from the second and third rows of Tab.1, the average decoded image quality for the 
five test images increased from 29.15 to 29.35 dB, whereas the proposed method provided 


better decoded image quality than Jacquin’s method. RBLVs contributed to 95% of the total 


13 


collage errors and played a key role in determining the total collage errors. However, 
compared with Jacquin’s method, the best-matched domain blocks of RBLVs in the proposed 
method were found in a larger domain block pool, EDBP, and smaller collage errors were 
obtained. The analyses are as follows: By Eqs. (4) and (5), for all range blocks, we define the 


following measure to compare the proposed with Jacquin’s methods: 


NumR NumR 
Measure=ACE?"P°s*4 = ACE Jacquin's — 3 CE (R,) = pa CE Jacquin's (R,) 24 
i=l 


i=l 
where ACE!ropesed and ACE*ain’s denote the ACEs of all range blocks for the proposed 
method and Jacquin’s method, respectively. If Measure is < 0, the proposed method can 
provide a smaller ACE and achieve a better block-matching effect than Jacquin’s method. 
Furthermore, for RBLVs, we define the following measure, Measure®®'Ys, to compare the 


proposed with Jacquin’s methods: 


M easure RBLVs _ ACE RBLVs Proposed _ ACE RBLVs_Jacquin's 


NumRBLV NumRBLV 25 
a9 RBLVs_Proposed RBLVs_Jacquin's 
= >) CE (R)- >) CE (R,) 


i=l i=l 


where ACERBLVs Proposed and ACERBLVs_lacauin's denote the ACEs of RBLVs for the proposed 
and Jacquin’s methods, respectively. NumRBLV denotes the number of RBLVs. Similarly, 
for RBSVs, we define the following measure, Measure®®SYs, to compare the proposed with 


Jacquin’s methods: 


Measure RBSVs =ACE RBSVs_Proposed _ ACE RBSVs_Jacquin's 
NumRB NumRBSV 26 


SV 
2 > CE RBSVs_Proposed (R,) = 5 CE RBSVs_Jacquin's (R) 


i 
i=l i=l 


where ACERBSVs Proposed and ACERBSVs_tacquin’s denote the ACEs of RBSVs for the proposed 


and Jacquin’s methods, respectively. NumRBSV denotes the number of range blocks with 


24-26 


large variances. From Eqs.** ~°, we have 


RBLVs RBSVs 27 


Measure=Measure +Measure 
For example, for the boat image in Fig. 7(a), we have Measure®®!Vs=- 14,359.43 < 0 which 
implies that for RBLVs, the proposed method can achieve smaller collage errors and provide 
a better block-matching effect. Although Measure®®SYS=13,751.68 > 0, by Eq., we have 
Measure=—607.75<0. Furthermore, by Eq.”4, we know that the proposed method can provide 


a smaller ACE, which implies a smaller ACER. Finally, by Eqs.* 5, we know that compared 


with Jacquin’s method, the proposed method can provide better-decoded image quality. 


(a) Boat (b) Peppers (c) Couple 


(d) Airplane (e) Lake 


Fig.8 Edge regions for the five test images. 
Further, in the proposed method, SDBP is mainly designed to improve the quality of 
RBLVs, which are mainly located in the edge regions. Moreover, compared with the smooth 
regions, the eyes of people are more sensitive to the edge regions, and the quality of the edge 


regions plays a crucial role in influencing people’s visual effects. Thus, we continue to 
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compare the quality of the edge regions of the previous methods with that of the proposed 
method. Figure 8 shows the edge regions for the five test images, and white pixels represent 
the edges. We used MSEE% in Eq. (28) to assess the quality of the edge regions in the 
decoded image regarding the counterpart in the original image, and a large MSE**2* implies 


poor image quality. 


1 NumE 


a 2) 
N i Edge __ Original Edge _ „„Decoded_Edge 28 
umE “7 


where NumE denotes the number of pixels in the edge region. pO“siral Edses and pPecoded_Fdgeg 
denote the pixel values of the same edge region in the original and decoded images, 
respectively. Tab.2 implies that for the five test images, the corresponding average MSE** 
decreases from 211.92 to 195.47. We know that compared with Jacquin’s method, the 
proposed method can provide smaller MSE “2°, which imply better-decoded image quality in 
the edge regions. 

From the fourth to seventh rows of Tab.1, we can see that compared with Jacquin’s 
method, the average percentage of the required computations can be reduced from 100 to 
79.72%, and the average bits per pixel of the proposed method can be reduced from 1.56 to 
1.32. Thus, we know that the proposed method can provide fewer bits per pixel and reduce 
the computations required. For Jacquin’s method, for an arbitrary range block, the 
corresponding block matching was restricted within conventional DBP. We had 
NumR=(M//B)x(N/B)=4096 range blocks and NumD’=(/2B)x(N/2B)x8=8192 domain 
blocks. In the encoding process, we should do block matching NumRxNumD’=33,554,432 


times. For the proposed method, the block-matching operations of RBLVs were conducted 


within EDBP. Because the EDBP consists of DBP and SDBP, which contain 
NumS=(//4B)x(N/4B)x8=2048 domain blocks, the EDBP contains 
NumED=NumD’+NumS=10240 domain blocks. For example, for the boat image in Tab.1, 
when BLAPACE reaches 95%, the NumRBLV=2728 range blocks are encoded, and we 
should do block matching NumRBLV xNumED=27,934,720 times. In the following encoding 
process, we adopted the no-search fractal encoding method for the remaining 1,368 RBSVs. 
Because the best-matched domain blocks of RBSVs are directly appointed, the corresponding 
block-matching operations can be saved, and the percentage of the computations required for 
the proposed method can be reduced to 27,934,720/33,554,432 = 83.25%. Contrarily, for 
each RBSV, we should only use 2 and 8 bits to quantize the scaling and offset coefficients, 
respectively. There is no need to store the isometric and position information of the 
best-matched domain blocks. Thus, the corresponding memory space can be saved. For 
example, for the boat image, the bits per pixel can be reduced from 1.56 to 1.35. 

Moreover, from Eq. (23), the proposed method can adaptively divide the range blocks 
into two categories, RBLVs and RBSVs. For different images, the numbers of RBLVs that 
contribute to 95% of the total collage errors vary significantly, and the corresponding bits per 
pixel and computations required fluctuate significantly as well. For example, the percentage 
of the required computations decreased from Peppers’ 90.97% to Airplane’s 55.48%, and the 
bits per pixel decreased from Peppers’ 1.42 to Airplane’s 1.13. 

Similar to Tables 1 and 2, Table 3 lists the performance comparison between 
Chaurasia’s and the proposed methods. The average decoded image quality increased from 


28.83 to 29.00 dB, average percentage of the required computations decreased from 100 to 
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64.89%, and average bits per pixel decreased from 1.56 to 1.33. Moreover, from Table 4, we 
observe that the average MSEF*2* decreases from 225.44 to 211.90. Clearly from Table 5, 
compared with Zheng’s method, the average decoded image quality of the proposed method 
increased from 28.96 to 29.16 dB, average percentage of the required computations decreased 
from 100 to 70.90%, and average bits per pixel decreased from 1.56 to 1.37. From Table 6, 
the average MSEF*2* decreased from 222.54 to 206.13. From Table 7, compared with Gupta’s 
method, the average decoded image quality of the proposed method increased from 29.03 to 
29.19 dB, percentage of the required computations decreased from 100 to 64.19%, and 
average bits per pixel decreased from 1.56 to 1.33. From Table 8, the average MSE ee 
decreased from 217.98 to 203.57. In summary, compared with the previous methods, the 
proposed method can provide a better quality of decoded images, require fewer bits per pixel, 


and reduce required computations. 


Table 1. Performance comparison between Jacquin’s method and the proposed method.!!! 


Test images Boat Peppers Couple Airplane Lake Average 
Jacquin’s 28.28 31.80 28.47 29.46 2712. 29.15 
PSNR(dB) i 
Jacquin’s+Proposed 28.39 32.03 28.60 29.65 28.06 29.35 
Computations Jacquin’s 100% 100% 100% 100% 100% 100% 
needed(%) Jacquin’st+Proposed 83.25% 90.97% 89.54% 55.48% 79.35% 79.72% 
Jacquin’s 1.56 1.56 1.56 1.56 1.56 1.56 


Bits per pixel : 
Jacquin’s+Proposed 1.35 1.42 1.40 1.13 1.32 1.32 


Table 2. Comparison of MSEF“®’s between Jacquin’s method and the proposed method."! 


Test images Boat Peppers Couple Airplane Lake Average 
Jacquin’s 250.32 123.04 204.69 217.08 264.45 211.92 
Jacquin’s+Proposed 235.47 111.97 191.17 203.41 235.31 195.47 


MSEE4¢ge 


Table 3. Performance comparison between Chaurasia’s method and the proposed method."! 


Test images Boat Peppers Couple Airplane Lake Average 


PSNR(dB) Chaurasia’s 28.00 31.47 28.01 29.37 27.28 28.83 


Chaurasia’st+Proposed 28.10 31.63 28.25 29.45 D757 29.00 


Computations Chaurasia’s 100% 100% 100% 100% 100% 100% 
needed(%) Chaurasia’st+Proposed 68.36% 75.15% 72.12% 44.09% 64.75% 64.89% 
Chaurasia’s 1.56 1.56 1.56 1.56 1.56 1.56 


Bits per pixel ; 
Chaurasia’s+Proposed 1.37 1.43 1.41 1.13 1.33 1.33 


Table 4. Comparison of MSE®’!°s between Chaurasia’s method and the proposed method."! 


Test images Boat Peppers Couple Airplane Lake Average 
Chaurasia’s 263.77 130.96 220.92 222.58 288.95 225.44 
Chaurasia’s+Proposed 251.42 123.25 206.89 212.98 264.96 211.90 


MSEEéee 


Table 5. Performance comparison between Zheng’s method and the proposed method.!! 


Test images Boat Peppers Couple Airplane Lake Average 
Zheng’s 28.04 31.78 28.00 29.29 27.70 28.96 
PSNR(dB) 
Zheng’s+Proposed 28.15 31.99 28.15 29.49 28.04 29.16 
Computations Zheng’s 100% 100% 100% 100% 100% 100% 
needed(%) Zheng’s+Proposed 85.74% 76.22% 77.48% 41.59% 67.46% 70.90% 
: . Zheng’s 1.56 1.56 1.56 1.56 1.56 1.56 
Bits per pixel 
Zheng’s+Proposed 1.52 1.41 1.42 1.14 1.36 ES i 


Table 6. Comparison of MSEF“’s between Zheng’s method and the proposed method.'! 


Test images Boat Peppers Couple Airplane Lake Average 
Zheng’s 265.78 123.69 230.40 226.85 266.00 222.54 
Zheng’s+Proposed 251.20 113.24 216.04 212.31 237.84 206.13 


MSEEéee 


Table 7. Performance comparison between Gupta’s method and the proposed method.” 


Test images Boat Peppers Couple Airplane Lake Average 
Gupta’s 28.16 31.75 28.30 29.29 27.66 29.03 
PSNR(dB) 
Gupta’s+Proposed 28.21 31.95 28.39 29.46 27.95 29.19 
Computations Gupta’s 100% 100% 100% 100% 100% 100% 
needed(%) Gupta’s+Proposed 67.38% 73.02% 73.00% 44.26% 63.28% 64.19% 
7 i Gupta’s 1.56 1.56 1.56 1.56 1.56 1.56 
Bits per pixel 
Gupta’s+Proposed 1.36 1.42 1.42 1.13 1.32 1.33 


Table 8. Comparison of MSEF“2¢s between Gupta’s method and the proposed method." 


Test images Boat Peppers Couple Airplane Lake Average 
Gupta’s 258.20 124.61 213.85 225.17 268.06 217.98 
Gupta’s+Proposed 247.13 114.19 202.64 213.12 240.75 203.57 


MSEFése 


5 Conclusion 


In this study, an adaptively hybrid fractal coding method was proposed to improve the 
performance of fractal coding methods. By analyzing the importance of the variances of 
range blocks, we found that RBLVs play a key role in determining the quality of decoded 
images. Thus, we designed EDBP for RBLVs to reduce the total collage errors. For the 
remaining RBSVs, the no-search fractal coding method, which can achieve fewer bits per 
pixel and computations required, was adopted. Moreover, if BLAPACE reaches 95%, an 
adaptive method for dividing the range blocks into two categories was proposed. In terms of 
the whole image, the proposed method improved the performance of fractal coding methods 


effectively and made them more useful in practical applications. 
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